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Abstract 

We construct families of cell complexes that generalize expander graphs. These families are 
called non-/c-hyperfinite, generalizing the idea of a non-hyperfinite (NH) family of graphs. Roughly 
speaking, such a complex has the property that one cannot remove a small fraction of points and be 
left with an object that looks k — 1-dimensional at large scales. We then consider certain quantum 
systems on these complexes. A future goal is to construct a family of Hamiltonians such that every 
low energy state has topological order as part of an attempt to prove the quantum PCP conjecture. 
This goal is approached by constructing a toric code Hamiltonian with the property that every low 
energy state without vertex defects has topological order, a property that would not hold for any 
local system in any lattice Z d or indeed on any 1-hyperfinite complex. Further, such NH complexes 
find application in quantum coding theory. The hypergraph product codesi of Tillich and Zemor 
are generalized using NH complexes. 
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In this paper, we construct several families of cell complexes which generalize expander 
graphs. More precisely, these families are a generalization of families of graphs which are 
"not hyperfinite" as defined by Elek— ; being not hyperfinite is a property that is closely 
related to being an expander graphs but slightly weaker. We then study various quantum 
Hamiltonians on these complexes and find that they realize properties that rely on the fact 
that these complexes are not k- hyperfinite, as we define below in definition 12.21 

The motivation for this work is that expander graphs play a key role in coding theory 
and complexity theory and we are looking for a similar construction which could play a 
similar role in quantum coding theory and quantum complexity theory. Later in the paper 
we discuss applications to coding theory, but now we briefly discuss the application to 
complexity theory, in particular the quantum PCP conjecture 14 . The quantum PCP (qPCP) 
conjecture deals with the difficulty of approximating the ground state energy of quantum 
many-body Hamiltonians. Suppose we have a system with N degrees of freedom, each degree 
of freedom having a Hilbert space of dimension 0(1) with the Hilbert space of the whole 
system being the tensor product of these Hilbert spaces. The Hamiltonian if is a sum of 
interaction terms, each interaction term being supported on 0(1) sites with each site having 
0(1) interaction terms supported on that site, and each interaction term having operator 
norm bounded by 0(1) so that the operator norm of H is 0(N). The qPCP conjecture is that 
there is some constant c (this constant may depend upon the particular quantities referred 
to as being "0(1)" in the previous sentence) such that it is QMA-hard to approximate the 
ground state energy to an accuracy better than cN. In Ref. Il2|, it was noted that in order 
for the qPCP conjecture to be true, there must exist a family of Hamiltonians H, obeying 
these requirements on the support and norm of the interaction terms, such that there is no 
trivial state ip (as defined below) such that (ip, Hip) is at most cN above the ground state 
energy. The reason that this needs to hold is that otherwise the trivial state can serve as a 
classical witness for the existence of a low energy state, putting the approximation problem 
wthin NP. Defining the "energy density" of a state to be 1/N multiplied by the difference 
between the energy of that state and the ground state energy, the requirement is that there 
exist no state with energy density less than c for some c. In this paper, we consider the easier 
problem of constructing a family of Hamiltonians with no low energy trivial states (NLTS). 
The NLTS problem, defined in definition II .3[ is a natural first step toward provng qPCP, and 
is in our view not a less interesting objective. The Hamiltonians that we construct in fact 
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will have a ground state energy that we can compute exactly, so this does not yet address 
more difficult problem of constructing a family of Hamiltonians such that approximating 
the ground state energy is QMA-hard. Our constructions use only Hamiltonians which are 
a sum of commuting terms. In fact, we do not fully succeed in constructing such a family 
of Hamiltonians but we obtain partial results which rely explicitly the complex being not 
1-hyperfinite. In particular we construct a complex on which a toric code has a certain "one- 
sided" NLTS, in that there is no trivial state with no vertex defects and with low energy for 
the plaquette terms, where throughout this paper we use "low energy" to mean an energy 
density bounded by some small constant. We conjecture a construction of a system with 
NLTS (removing this "one-sided" restriction) but we do not prove this. 

This application to quantum systems led to our definition 12.21 which is different from 
other generalizations of expander graphs considered elsewhere^ 1 ^. The outline of the paper 
is as follows. In section [1] we define trivial states and quantum circuits. In section [21 we give 
the definition of fc-hyperfiniteness and for every k we construct families of complexes which 
are not /c-hyperfinite. In section [3], we further expand on the motivation for considering these 
complexes, by giving an analogy to classical systems. Sections l2|3l can be read in either order; 
more mathematical readers might prefer our order, while more physical readers might prefer 
to read section [3] first. Section [3] can be skipped if desired as it is purely motivational, 
however it does present some results about statistical mechanics on expander graphs that 
may be of interest; while these results in this section may be well-known, we have not found 
them in the literature. In section H] we review the definition of the two-dimensional toric code 
Hamiltonian, and give this definition for arbitrary 2-complexes. Further, we describe other 
toric codes, so-called (q, r)-codes in dimension d = q + r (often, these would be referred to as 
(p, q) codes, but we are using the symbol p for other purposes here). In section [5l we consider 
the toric code on a family of complexes that is not 1-hyperfinite. While we do not show that 
the code on this family has NLTS, we show that there is no trivial state for which the terms 
in the Hamiltonian that penalize broken strings are in their ground state and which is low 
energy with respect to the other terms (the terms that give a dynamics to the string). As 
we note, this property would not occur for any complex which was 1-hyperfinite including 
such familiar examples as a two-dimensional lattice or an expander graph. In section [6] 
we discuss the relation to hypergraph product codes and give an alternative presentation 
of those codes in terms of the language of chain complexes; this presentation leads us to 
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a generalization of these codes. We also present a possible conjecture on a Hamiltonian 
with NLTS. Finally, in the discussion we discuss an interpretation of our results in terms of 
dynamics of strings or higher dimensional surfaces, and give a conjecture that if a certain 
family of non-/c-hyperfinite complexes could be constructed then one could indeed give a 
family of complexes on which a toric code hadNLTS, now requiring only that the state 
be low energy with respect to all terms in the Hamiltonian, rather than requiring that it 
also be zero energy with respect to certain terms. In an appendix, we discuss possible 
generalizations of our definition of hyperfiniteness to the case of noncompact spaces where 
we define a notion of A;-amenability. 



1. DEFINITIONS AND TRIVIAL STATES 

We define some metric dist(i, j) between degrees of freedom We choose this metric so 
that the interaction terms in the Hamiltonian have support on sets of diameter 0(1). One 
way to define this metric would be to define an interaction graph, with vertices corresponding 
to degrees of freedom and an edge connecting two vertices if both degrees of freedom were 
in the support of some term in the Hamiltonian, and then to use the metric which is the 
shortest path metric on this graph. In fact, later in the paper we will not use this metric 
but a different one; the reason is that for much of the paper, the degrees of freedom will be 
on edges, rather than on vertices. So, for that situation we will choose a metric in which two 
edges attached to the same vertex are a distance 1 from each other. However, the specific 
metric used is not of much importance (the mathematical notion of quasi-isometry expresses 
in detail which metrical dfferences are "unimportant" for us). 

We begin by defining a quantum circuit. 

Definition 1.1. A unitary quantum circuit U is a unitary of the form 

U = U 1 U 2 ....U k , (1) 

where k is called the "depth" of the circuit, such that each unitary Ui can be written as a 
product 

Ui = U itl Ui,2...U itrH (2) 

where the support of U i)CL is disjoint from the support of for a ^ b for all i and where 
the integer is the number of unitaries in the i-th round. The "range" of a unitary U i>a is 
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the diameter of the support of that unitary. We define the range of a quantum circuit to be 
the product of the depth of the quantum circuit times the maximum range of the unitaries in 
that quantum circuit. 

Note that this definition of the range of a quantum circuit is not completely standardized 
since some authors take the range of the circuit to be the maximum range of a unitary in 
the circuit; we prefer this definition here since it makes the definition of trivial states more 
natural. 

We define a trivial state as follows: 

Definition 1.2. A state ip with \ip\ = 1 is (R,e)-trivial if there exists a unitary quantum 
circuit U with range R such that 

\ip - Utp prod \ < e (3) 

for some product state ip pro d- 

We say that ip is R-trivial if it is (R, 0) -trivial. 

One may allow the use of ancillas in these constructions. For each degree of freedom in 
the original system, we define two Hilbert spaces, one called the real Hilbert space which 
is the same as the Hilbert space for that degree of freedom in the original system and one 
called the ancilla Hilbert space. We can then consider trivial states on this larger Hilbert 
space obtained by unitary quantum circuits acting on product states on the larger Hilbert 
space. The ancillas are organized by the same geometry as the original Hilbert space; that 
is, the distance dist(...) between an ancilla on a given site and the real degree of freedom on 
that site is equal to zero. By tracing out the ancilla degrees of freedom, the result may be a 
mixed state; we refer to this as a trivial mixed state, and indeed any convex combination of 
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Our results later are 



such states will be referred to as a trivial mixed state following Ref. 
phrased in terms of trivial pure states, but the same results can be proven for trivial mixed 
states (add ancillas to prove it for a trivial mixed state obtained by tracing out ancillas and 
then once we have proven that no such trivial state has low energy it will also follow that 
no convex combination of such states will have low energy). 

With these definitions, we are now prepared to define the NLTS problem. 

Definition 1.3. Consider a family of Hamiltonians, H^, with N — > oo. Each acts on 
a Hilbert space which is a tensor product of N different finite dimensional Hilbert spaces, 



called sites, each having dimension 0(1). Each is a sum of interaction terms, each 
interaction term being supported on 0(1) sites with each site having 0(1) interaction terms 
supported on that site, and each interaction term having operator norm bounded by 0(1). 
Assume that the smallest eigenvalue of is equal to for all N . 

We say that such a family of Hamiltonians has NLTS if there exists a constant e > 
such that there is no R such that for all N there is an R-trivial state ipN such that 

(if> N , H N ifj N ) < eN. (4) 

We conjecture 

Conjecture 1.4. A family of Hamiltonians with NLTS does exist. 

We have not specified that Hn should be a sum of commuting terms in the above defini- 
tion, or that H]\f should be frustration free, although the Hamiltonians that we consider n 
the rest of the paper are indeed frustration-free and sums of commuting terms. One might 
naturally defining problems such as commuting NLTS where such a commuting restriction 
is imposed. 

Our results later will show the absence of .R-trivial states with certain low energy prop- 
erties. However, from the point of view of the qPCP conjecture and NLTS problem, these 
definitions of (R, 0)-trivial and (R, e)-trivial are not too different. If there is an (R, e)-trivial 
state if), then 

\(Uilj prod ,HUilj prod } - (ip,Hif;}\ < (2e + e 2 )||#||, (5) 

and so, since ||fZ"||/./V = 0(1), if there is no .R-trivial state for a Hamiltonian with energy 
density less than c, then there is no (R, e) trivial state for that Hamiltonian with energy 
density less than c + 0(l)(2e + e 2 ). This point may not apply to our one-sided results later, 
since they are based on states with zero vertex defects (hence zero, rather than small, energy 
for a certain type of term), but still we find the case of .R-trivial states worth considering. 

We note that the reason that trivial states are relevant for the qPCP conjecture is that a 
trivial state can serve as a witness because one can efficiently compute the reduced density 
matrix of the trivial state on small sets. However, NLTS is weaker than qPCP because there 
are lots of other possible choices for a classical witness. For example, one could take a state 
that is trivial using a different metric to define the support of the unitaries Ui )a in Eq. (J2J), 
different from that used to define the interactions. We will not consider such states in this 
paper. 
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2. FAMILIES OF COMPLEXES THAT ARE NOT fc-HYPERFINITE 



In this paper we will construct various complexes. In Ref. [12|, the complexes constructed 
were simplicial complexes. In this paper, it is also convenient to use cubical and more 
generally polyhedral complexes where convex polytopes are glued together via simplicial 
homeomorphisms of their faces. Another distinction is that in Ref. [12I a complex called an 
interaction complex was constructed: the degrees of freedom of the quantum Hamiltonian 
were placed on the 0-cells of the complex, and the higher cells of the complex were determined 
by the support of the terms in the Hamiltonian. In this paper, we will often place the degrees 
of freedom on 1-cells or even high cells of the complex. 

A graph can be viewed as a 1-complex. We will use the terms "vertex" and "0-cell" 
interchangeably; we will similarly use the terms "edge" and "1-cell" interchangeably and 
also the terms "plaquette" and "2-cell" interchangeably. We will later refer to "vertex 
operators" and "plaquette operators" for a toric code, which will be operators defined for 
each 0-cell or 2-cell of the complex, respectively. One case in which we will try to distinguish 
between these terms is that for certain complexes which are constructed from a product of 
graphs, we will use the term "1-cells" for the complex and the term "edge" for the graph. 

We place a metric on a 1-complex, by defining each edge to have length 1, and using a 
path metric. This reproduces the usual graph metric between vertices on the graph (the 
distance between neighboring vertices is equal to 1). We extend this metric to some path 
metric over the higher cells so that the cells in each dimension have bounded diameter. Now 
define fc-localizable complexes by: 

Definition 2.1. A metrized simplicial l-complex Ki is "k-localizable with range R" if there 
exists a continuous function f from K\ to some simplicial k complex Kk (using the same 
graph metric as above on the edges of Kk, extended as above to the higher cells of Kk) such 
that the diameter of the pre-image of any point in K k is less than R. 

Next define /c-hyperfiniteness: 

Definition 2.2. Consider a family of I- complexes, Ki(N), where N is the number of 0-cells 
in Ki(N), with N — > 00. Such a family is said to be "k-hyperfinite" if for all e > 0, there 
exists an R such that for all N one can remove at most a fraction e of the 0-cells of Ki(N), 
while removing also all attached higher cells, such that the resulting complex is k-localizable 
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with range R. 

Families of interest to us generally have bounded local geometry: the number of j ' + 1 
cells meeting any j-cell and the number of j cells in the boundary of any j + 1-cell is upper 
bounded independent of N (and j). But we do not include this property in the definition. 

The definition of a family of 1-complexes being O-hyperfinite coincides with the definition 
of a family of graphs being hyperfinite^, motivating our terminology. 

Note that if the metrics on the family Ki(N) are altered by a uniformly bounded 
amount: d N replaced with a quasi-isometric d' N so that for all x, y e Ki(N), -^djy(x,y) — 
C\ < d' N (x,y) < C d N (x,y) + C\ for C , G\ > 0, the property of being or not being k- 
hyperfinite is unchanged, so the choice of metric is quite flexible. 

We now construct families of complexes which are not Z-hyperfinite. To exposite key ideas 
first we begin with the case 1 = 1. Consider a family of expander graphs with fixed degree d 
and with girth g of the graph diverging as iV — > oo; in fact, we may take the girth of order 
log(iV). A graph G in this family with N vertices has edges and has first Betti number 
bi(G) = ^piV+l. We will define a 2-complex, by taking the Cartesian product of this graph 
G with itself, G x G. This is the "product graph" with 2-cells, i.e. Euclidean squares, of 
the form edge x edge glued in. Note that a 0-cell in G x G is labeled by a pair (i,j), where 
i,j are vertices in G. A 1-cell in G x G is labeled by either a pair (e, i), where e is an edge 
in G and i is a vertex, or a pair (i, e). Since the girth of the graph G is diverging, we may 
assume that it is larger than 4 for large enough N, so all length 4 loops are boundaries of 
such squares. Explicitly, let an edge e connect vertices e%, t2 in G and let an edge / connect 
vertices fx, f2- Then, the length 4 loops are of the form (e, fx), (ei, /), (e, fa), (f, e^). An 
important property of G x G is that it has N 2 0-cells and has an "extensive" second Betti 
number, where a quantity will be called "extensive" if it is 0(the number of 0-cells). In this 
case, 



It is this large second Betti number, combined with the high girth of G, that we use to prove 
that the family of complexes G x G is not 1-hyperfinite. 

Before giving the detailed proof of the theorem below, let us sketch the overall strategy. 
First, we use the fact that the second Betti number of G x G is extensive to show that 
deleting a sufficiently small fraction of cells from G x G leaves a complex G x G \ C that 




(6) 
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still has a nonzero second Betti number since deleting a 2-cell changes the Betti number by 
at most one. This implies that the resulting complex still contains a nontrivial 2-cycle. The 
assumption that G x G \ C is 1-localizable means that there is a map / from G x G \C 
to a 1-complex K. This map / maps the nontrivial 2-cycle to a subset of K. It is at this 
point that we use the large girth assumption on G to show a contradiction: since G x G 
has large girth, the nontrivial 2-cycle in G x G \ C cannot be continuously mapped to a 
1-complex using a map with bounded diameter of pre-images. To motivate why the large 
girth is needed, consider the following case of a complex with a small girth: a "long, thin" 
torus (a torus with the second radius being small; this torus is the product of a large girth 
line graph with a small girth line graph). One can find a continuous map from this torus 
to a circle in the natural way, mapping a point on the torus with angles (9, <fi) to a point 
on the circle at angle 9, giving a map with a diameter of pre-image bounded by the girth of 
the second graph. In contrast, since we are considering complexes G x G with large girth, 
we will show that a continuous map to a 1-complex with bounded diameter of pre-images 
is impossible, showing that such a map would lead to a contradiction as it would instead 
imply that G x G\C had vanishing second Betti number. 

A second point worth recalling is the idea of a simplicial map and the simplicial approx- 
imation theorem. These ideas enable us to avoid pathologies that might occur for arbitrary 
continuous functions and let us instead deal with more combinatoric ideas. Simplicial maps 
are maps such that the images of the vertices of a simplex span a simplex. Simplicial maps 
are completely determined by their action on the 0-cells. The simplicial approximation the- 
orem enables us to approximate a continuous map between two complexes by a simplicial 
map, after some subdivision of the simplices of both complexes. The approximation can be 
made arbitrarily accurate by taking sufficiently fine subdivision and the approximating map 
is homotopic to the original map. 

Theorem 2.3. Let Gjv be a family of connected degree = d graphs with diverging girth. The 
family of 2- complexes G^ x G^ is not 1-hyperfinite. 

Proof. To a connected graph G is associated a Z 2 -cellular chain complex: 

Ci A C A Z 2 -)• 

where Cj is the Z 2 - vector space generated by i-cells: C\ = Z 2 2 , C = Z£ , and the final 
map (epimorphism) to Z 2 is the "augmentation" a which records connectivity of the graph 
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by taking every generator to 1. The map d is onto ker(a), so 61 = rank(i?i(G; Z 2 )) = 
rank ker(<9) = ^ + 1. By the Kiinneth formula, b 2 (G x G) = rank H 2 (G x G; Z 2 ) = 
(^N+l) 2 . 

Deleting a 0-cell implies deleting d 2 adjacent squares (2-cells) of G x G. Since d > 3 is 
fixed, for iV large and e < 0.02 which we now assume, deleting an e-fraction of 0-cells implies 
deleting at most d 2 eN 2 < {^N + l) 2 2-cells from G x G. 

An easy exercise in Mayer- Vietoris (M-V) sequences establishes: 

Lemma 2.4. b 2 {G x G \ C) > {^N + l) 2 - \C\ > (^N + l) 2 - d 2 eN 2 > 0, where 
C = C\ U C 2 U ■ ■ • U C\c\ is the union of the 2-cells being deleted and \C\ is the number of 
such 2-cells. Clearly for e small enough b 2 {G x G \ C) becomes extensive. 

Proof. If A and B are sub-complexes of a complex X the M-V exact sequence reads 

— > Hi(A HB)^ Hi(A) © H t (B) — > H t (X) UB)^ 

where i is arbitrary and we take all coefficients = Z 2 . Let us set i = 2 and inductively 
assume I = GxG\(7 1 U---U C fc _i, A = C fe , and B = G x G \ Ci U • • • U C k . Now the right 
hand term H\(A D -B) is the first homology of some subset of the circle dCk which is Z 2 if 
the subset is the entire circle and otherwise. H 2 (Ck] Z 2 ) = 0, a disk being contractible has 
no homology except in dimension 0. So our sequence reads: 

— )• H 2 (G x G \ Ci U • • • U C k - Z 2 ) — )• B 2 (G x G \ Gi U • • • U C fc _i; Z 2 ) — )• Z 2 or — )■ 

By exactness, the rank of the left term, b 2 (G x G \ Gi U • • • U G*,; Z 2 ), plus the rank 
of the right term must be greater than or equal rank(B 2 (G x G \ C\ U ■ ■ • U G^-i; Z 2 )) = 
6 2 (G x G \ Gi U • • • U Gfe_i; Z 2 ). Thus the Z 2 second Betti number b 2 will either decrease by 
one or remain unchanged whenever a 2-cell C k is deleted. The lemma follows. □ 

Lemma 2.5. The map H 2 (G x G \ G; Z 2 ) H 2 (G x G; Z 2 ) an injection. 

Proof. The map inc 2 fits into the exact sequence of a pair: 
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— ► H 3 (G x G,G x G\C;Z 2 ) # 2 (G xG\C;Z 2 )^ # 2 (G x G; Z 2 ) — > 

||? excision 
tf 3 (C,<9C;Z 2 ) 

Ik 


Via excision, the first term is H3 of a relative 2-complex and therefore vanishes, showing 
inc 2 to be an injection. □ 

Remark 2.6. The vanishing of homology in dimensions exceeding the dimension of the 
complex is a property of all ordinary homology theories, i.e. theories obeying the dimension 
axiom. These include H*( ; Z 2 ) which we are using but not, for example, theories such as 
K -homology, ). 

Proof of theorem. Pick e > small enough so that for N large and for any set C of 2-cells 
which may be deleted the Z 2 -Betti number b 2 (G x G \ C) is extensive (actually, non-zero 
suffices for the proof). Write X := G x G\C. Suppose X is 1-localized with range R (for all 
large N) and that this property is witnessed by a map / : X — )■ K, a 1-complex, i.e. a finite 
graph. It is convenient (but not at all essential, as we will see in the proof of Theorem 12 . 71) 
to represent a non-zero 2-cycle in H 2 (X; Z 2 ) by a map g from a (not necessarily orientable) 
surface S. 

g : £ — > X satisfies: 

<7*[£] 2 = x ^ E H 2 {X- Z 2 ) 

Now consider the pre-images W p '■— (f o g)^ 1 (p), p G K, of the composition. We know 
that for non-empty W p , diam(g(W p )) < R. At this point to avoid pathologies we will assume 
(permitting possible subdivision of the cell structure of all spaces E, X, and K) that / and 
g are simplicial maps. If a continuous map h between compact metric spaces has pre-images 
of diameter < r then Ve>03<5>0 such that V points p, diam(/?, _1 (Ball e (p))) < r + 5. 
The simplicial approximation / to / can satisfy ||/ — /|| SU p < e for any fixed e > 0. Since 
/ _1 (p) < f ~ l (,B e {p))i we have diam(/ _1 (£>)) < r + 5. Consequently, restricting to simplicial 
maps / has no effect, even for the constant, on the property of being "i?-localized" . 
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Let us give the idea of the contradiction we are approaching. The surface £ is "foliated" 
by "not-very-large-in-X" sets W p , diam g(W p ) < R. But 



)IcGxG 
* 



m 

cs 



and the Cartesian geometry d = \Jd\ + d\ on G x G built from the 1-complex path metric i 
each factor is locally convex, meaning it has convex neighborhoods — whose radius approach 
infinity as iV — > oo — while R is fixed. This suggests that the map W p — >■ G x G extends 
over the mapping cylinder M: 



M 



■GxG 



where M = £ x [0, l]/o~i x 1 = lt 2 x 1 whenever / o g{a{) — f o g(a 2 ). The rough thought 
is that M is "foliated" by cones to various points in K := image {fog) C K and that 
each of these cones separately maps (even canonically to its centroid) into GxG using the 
convexity structure. So far so good, but we would like to state that these cones fit together 
continuously so that h above is a continuous map. If we had a continuous h we would have 
a contradiction, since the induced diagram (below) on H 2 ( ; Z 2 ) shows that x — 0, whereas 
by hypothesis i ^ G H 2 (X; Z 2 ). 




inc* 

H 2 {M-Z 2 ) 



x^OE H 2 (X; Z 2 ) 

h* 



a.= inc* 



y^0EH 2 (GxG;Z 2 ) 



H 2 (K ;Z 2 )=0 

Lemma [2751 and i ^ imply y := a(x) ^ 0, but, in contradiction, the top map factors 
through H 2 (M; Z 2 ) = 0. To see this vanishing, note that any mapping cylinder deformation 
retracts to its target. In our case, this means the inclusion i : Kq C M, induced by 
Exl c Ex[0,l], has a homotopy inverse n induced by the projection p : £ x [0, 1] — > E X 1; 
ix oi ~ idx , and ion ~ id^v/- Consequently, and M are homotopy equivalent, K — M. 
By Remark l2.6[ since Kq is a 1-complex, H 2 (Kq; Z 2 ) = 0, a property inherited by H 2 (M; Z 2 ), 
since homotopy equivalences induce isomorphisms on homology. 
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The preceding outline is almost complete. The flaw is that pre-images of points even under 
a continuous map can jump discontinuously (consider a not-strictly monotone function from 
R to R which is constant on exactly one closed interval. All point pre-images are single 
points except for that one interval.) So our proposed "cone to centroid" construction might 
produce discontinuities in h. But these are easily repaired. 

It was already noted that the pre-images under / of e-balls also satisfy diam f^ 1 (B e (p)) < 
R for sufficiently small e > 0. This allows us to subdivide the target graph Kq so finely 
that the pre-image of every 1-simplex (as well as every 0-simplex) has diameter < R. Now 
build a slightly different (in topology it would be called the "block" version) but homotopy 
equivalent version M' of the mapping cylinder M and build a map from M' to G x G by 
using the local convex structure of G x G twice. First cone (/ o g)~ l (0-cells) and map these 
cones linearly to the centroid of g(f o g)^ 1 (0-cell) C G x G. Then for each 1-cell e C Kq 
cone: (/ o g)~ x {e) U cone((/ o g)~ 1 (d_e)) U cone(/ o g)~ l (d + e) (Note: the latter two terms 
of the union and their maps to G x G were just previously constructed.) and again map 
this "secondary" cone to G x G linearly to its centroid, M' is easily seen to be homotopy 
equivalent to M. The easiest way is to construct a deformation retraction M' — > Kq: first 
compressing the fibers above 0-cells of K and then using Urysohn's Lemma to deform the 
second cone in a neighborhood of (f °g)~ 1 (de) to a mapping cylinder structure as illustrated 
in Figure 2.1 below. 




t 



v 



original cone of cones 



deform map here deform map here 

final mapping 
cyclindary structure 



FIG. 2.1. Rotate about x-axis to make illustration 3D 



This completes the proof of Theorem 12. 31 □ 

There are immediate generalizations of Theoren l2.3l in several directions. One direction 
is to find other families of spaces, e.g. manifolds with the local geometry of a symmetric 
space which, for topological reasons, have extensive Z 2 homology in some dimension k, 
— — > constant, and to prove that these families are not (k — l)-hyperfinite. But the most 

volume — c \ j j sr 
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straightforward generalization, the one presented next, is to simply consider fc-fold product 
of large girth graphs. 

Theorem 2.7. Let be a family of connected degree = d graphs with diverging girth. The 
family of k complexes G k N := Gn x • • • x Gn (k-copies) is not (k — 1 )-hyperfinite. 

Proof. The argument is parallel to the proof of Theorem 12.31 so we only highlight the differ- 
ences. 

Whereas the second Z 2 -Betti number b 2 was extensive for GxG, an iterated application of 
the Kiinneth formula shows that the (only) extensive Betti number for G k is = (^-N+l) k . 
A very similar homological argument shows that for e > small enough deleting an "e- 
fraction" of G k results in a subspace X whose b^ is still extensive. 

Once k > 2 an essential cycle x G H^{X] Z 2 ) may no longer be represented by a map of a 
/c-manifold, so our surface £ must be generalized to a fc-complex S fc with iffc(S fc ; Z 2 ) = Z 2 
generated by so that the class ^ x = g*[^ k ] for some g : S fc — > X. The fact that 
S was a surface, not a general 2-complex, actually played no role (other than an aid to 
visualization) in the proof of Theorem I2.3[ so we are not now inconvenienced by needing to 
take E fc to be a complex. 

Again for contradiction we assume there is a witness to (k — l)-hyperfmiteness in the 
form of a map / : X — > K, to a (k — l)-complex. Again we think of E fc as "foliated" by 
pre-images (/ o g)^ 1 (p), p G K, and diam(g(/ o g)^ 1 (p)) < R. G k still has a local convex 
Cartesian geometry so we should attempt to factor the map inc o g : £ fc — y G k through M, 
the mapping cylinder. S fc x [0, 1]/ (a, 1) = (a', 1) iff f(a) = f(cr'), as follows: 

E fc 9 - X inc - G k 

h ^ ' " 

\l 

Kq C K, a k — 1 complex 

All the homological facts are precisely parallel replacing H 2 ( ; Z 2 ) by ; Z 2 ), specifically 
since image (/) := K C K is a (k— l)-complex so Hk(M] Z 2 ) = 0, again giving a homological 
contradiction to the existence of k th homology in X. 

In higher dimensions (k > 2) a bit more needs to be said about the construction of M', the 
"block" model for M. Now instead of two stages of coning there will be k successive stages. 
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Triangulate Kq finely so that diam(/~ 1 (any cell) ) < R. Consider first each / _1 (0-cell) and 
attach a cone on each of these to S fc . The map h ■— incog : £ fe — > G k should be extended over 
these cones by "cone to centroid" using the local convexity structure of G k . Next consider 
each / _1 (e), e a 1-cell, and attach a cone to (/ _1 (e) U cone f~ 1 (d_e) U cone / _1 (<9 + e)). 
Again "cone to centroid" extends h over this secondary cone. Then consider each f -1 ^), 
C2 a 2-cell, and attach a cone to: 

r 1 (e 2 )U( |J coneCr 1 ^) U cone f-\d_e x ) U cone r\d + e x ))) 

eiG9e 2 

and extend h by "cone to centroid". Continue in this way until each cell of the final 
stratum is assigned a cone: cone(/ _1 (efc) U (lower cones)), and that cone is mapped via 
"cone to centroid" into G k . Again an elementary "bending of cone lines" shows that this 
block version M' is also homotopy equivalent to K and therefore has vanishing H^(M'; Z 2 ). 
This contradiction completes the proof. □ 

There are numerous reasons detailed below related to efficiency of quantum codes and to 
qPCP not to be satisfied with product of graph examples but to exploit the full power of Lie 
theory, Riemannian geometry, and topology to study more exotic contexts in which non-(/c — 
l)-hyperfinite manifolds arise. The duality reflected to the toric code fits very clearly with 
Poincare duality in manifolds and makes the search for non-(fc — l)-hyperfinite manifolds, 
rather than merely complexes, important. Surprisingly, we find that the mathematical 
questions that attempts at such constructions raise are just beginning to be addressed in 
the geometry literature. In particular, the calculation of L 2 -cohomology for symmetric 
spaces has surprising relevance. 

Let us take a brief tour of some intriguing classes of examples and record what is known, 
and what is open for these. 

Just as the hyperbolic plane H 2 is in some ways a homogeneous analog of an infinite 
(i-valent tree, with d being the analog of the Gaussian curvature, closed hyperbolic surfaces 
with injectivity radius r approaching infinity (meaning every point p G S enjoys an embedded 
ball neighborhood B r (p) C E of radius r) are analogous to a family of <i-valent graphs with 
increasing girth. (And further, surfaces such as those of arithmetic type, with Ai of the 
Laplacian bounded away from zero, are analogs of expander graphs.) By arguments utterly 
parallel to Theorems 12.31 and 12.71 we may prove: 
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Theorem 2.8. Let be a sequence of compact hyperbolic surfaces with injectivity radius 
r going to infinity as a function of the area A (r w | log e A is very nearly the optimum). 
Then the k-fold Cartesian products form a non-(k — 1) -hyper finite family. 

Proof. Apply the earlier arguments to the k th homology. is extensive. In 

fact, combining the Gauss-Bonnet theorem and the Kiinneth formula one computes that 
lim^oo &fc(S fc )/vol 2 fc(S fc ) = jp. A technical point is that our earlier discussion calibrated 
Betti number in terms of the number of cells in a complex and now we calibrate in terms of 
volume. Fortunately the following theorem of Ref. i3J allows a direct translation. 

Theorem 2.9. Given an integer d > 2 and a real number r > 0, there is constant C(d, r) > 
so that every hyperbolic d-manifold with injectivity radius > r can be triangulated with 
geodesic d-simplicies o~i of bounded geometry in the sense that each admits a homeo- 
morphism hi : o~i — > o~ , er the hyperbolic simplex with all sides of length 1 so that 

C [a, r) 

for all x,y G a; t . 

Proof. See Ref. S □ 

Given this theorem, the modification from the proof of Theorems 12.31 and 12.71 is routine. 

□ 

According to the Chern-Weil theory, the Euler characteristic x of a Riemannian manifold 
may be computed by integrating the Pfafnan of the Levi-Civita curvature over the manifold. 
For closed odd-dimensional manifolds the result is zero and for hyperbolic 2n-manifolds M 2n 
a local expression shows that there are dimension-dependent constants c n so that the Euler 
characteristic satisfies: 

X (M) = c n vo\(M). 

ci = — ^ and in general the c dd are negative and the c even positive 5 . 

Since there is also a combinatorial formula for x(M) in terms of the real (equivalently 
rational, Q) Betti numbers: 

X (M) = nU-l) l bl 

We learn that certain sums must be extensive. For example, we have this table for closed 
hyperbolic 2n-manifolds: 
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dimension = 2n 


extensive quantities 


2 


bf 


4 




6 


6? + b® + b®, bf + bf 


8 


bf + bf + bf, bf + bf 



(Poincare duality implies fej: = bf n _ k , accounting for the multiple extensive quantities 
listed.) 

In all our application either 0(l)-depth quantum circuits or length scales protecting 
quantum information arise so we are really only interested in toric codes on sequences of 
manifolds whose injectivity radius diverges. In this case there is more refined information^ 
coming from the behavior of L 2 -cohomology under the processes of taking geometric limits 
of symmetric spaces. 

Theorem 2.10. If {M d } is a sequence of hyperbolic d-manifolds, d fixed, with injectivity 
radius diverging, then all densities bf / vol — > except bf / vol, which approaches positive 
universal values for d = 2n. 

In building quantum codes Z p Betti numbers are more relevant than rational ones; the 
toric code is based, for example, on Z 2 -homology. Fortunately the universal coefficient exact 
sequence for arbitrary X: 

Hi{X- Z)®Z P -^ Hi(X- Z p ) — ► Tor(^!(X; Z), Z p ) — ► 

II* 

torsion(i7j_iX; Z) <S> Z p 

together with the identity 

rank(Free(^(X; Z))) = bf(X) 

shows that the mod p Betti numbers satisfy: 

bf p (X) = rank(#;(X; Z p )) > rank(Free(^(X; Z))) = bf{X). 

So bf serves as a lower bound for all bf p . 

Again a straightforward generalization of the proofs of Theorems 12.31 and 12.71 establishes 
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Theorem 2.11. Any sequence {M} of hyperbolic 2n-manifolds with injectivity radius di- 
verging is non-(n — 1 )-hyperfinite. 

Note that all complex matrix groups are residually finite, and by Selberg's lemma each 
admits torsion-free subgroups of finite index. These two facts permit the construction of a 
sequence of hyperbolic manifolds with diverging injectivity radius within the covering space 
tower of any fixed hyperbolic manifold. 

Although we now have excellent control of rational Betti numbers we really only have 
lower bounds on Z p Betti numbers. We have an: 

Open Question 2.12. For {M} a sequence of hyperbolic d-manifolds, d fixed = 3, 4, 5, 6, . . 

with injectivity radius diverging, is it possible that a Betti number bf p , for i ^ ~ could have 
nonzero density in the limit? 

Discussion: This question is number theoretic in nature. Investigation of analytic pro-p 

z I 

groups suggests some scaling like 6 1 p ~ vol 2 might be achieved for certain sequences at 
arithmetic 3-manifolds (with diverging injection radius). 

Curiously convincing heuristics and copious numerical evidence^ support the conjecture 
that for certain sequences of arithmetic hyperbolic 3-manifolds (with diverging injection 
radius) , 

log (order torsional (M; Z))/vol — > — , 
being the L 2 -analytic torsion of H 3 . 

Unfortunately, the numerics show that many primes and high prime powers factor this 
order so this observation does not imply extensive b l p for any p. See the discussion for the 
relevance of this open question. 

3. QUANTUM AS THE SQUARE OF CLASSICAL 

Before launching into the definitions and constructions, we give some additional motiva- 
tion for this problem. As noted, for the qPCP conjecture to be true, we need a family of 
Hamiltonians with NLTS. Sometimes, the property that a state is not trivial is referred to 
as a state being "topologically ordered" , although often the term "topological order" is used 
in a slightly different sense. Let us return to the classical setting and consider the property 
of long-range order: a long-range correlation between two degrees of freedom in the ground 
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state. We begin with one-dimensional systems with long-range order in the ground state, 
but not at thermal equilibrium in non-zero temperature. Then, we turn to two-dimensional 
systems, with long-range order at sufficiently small non-zero temperature, but which have 
low energy states without long-range order, and finally we discuss these systems on expander 
graphs, showing that every low energy density state has long-range order in a certain sense 
(this property may be known elsewhere but we do not know a reference to it in the litera- 
ture). We then give a parallel discussion in the quantum case; the graphs that we consider 
in the quantum case are precisely the square of the graphs we considered in the classical 
case. We note that the two-dimensional toric code has a nontrivial ground state, but has 
trivial states at nonzero temperature. We then turn to a four dimensional-toric code which 
has topological order for sufficiently small non-zero temperatures but which has low energy 
states which are trivial. This then motivates a consideration of a toric code on the square 
of an expander graph, which we consider in the next section, section [5j 

Consider the properties of the classical ferromagnetic Ising model on various graphs. To 
define this Ising model, for each graph G, we have a spin-1/2 degree of freedom on every 
vertex. The Hamiltonian is 

Hp erro = — ^ ] St Sj , (7) 

<i,j> 

where 5Z<» j> denotes a sum over all pairs of vertices i,j which are neighbors on the graph. 
If the graph is a one-dimensional line, then the Ising model has long-range order in the 
ground state: while there is a two-dimensional space of ground states, either all spins up or 
all spins down, the expectation value (if), S-Sjip) is equal to plus 1 in all ground states if), for 
all sites regardless of whether or not i, j are neighbors. However, the long-range order 
in this case is not stable to non-zero temperature. Let p(T) = exp(— i?/T)/tr(exp(— H/T)), 
so that p is the thermal density matrix at any temperature T > 0. The expectation value 
tr(p(T)S?Sj) decays exponential in the distance between i and j for any T > 0. However, 
if we move to a two-dimensional Ising model, the long-range order does survive at small 
enough non-zero temperature. Consider a families of graphs Gl, where Gl is a square- 
lattice of size L-by-L, and let Pl{T) be the corresponding thermal density matrix for graph 
Gl- Then, there is a T c > such that for all T < T c there is a quantity m(T) > such that 
tr(pL(T)SfSj) > m{T) for all L and all vertices i,j in Gl- 

However, while the long-range order survives at non-zero temperature in the two- 
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dimensional (and higher) model, there are states ipi ow of arbitrary low energy density which 
are not long-range ordered. We define the energy density of a state as follows: 

Definition 3.1. Consider a Hamiltonian H written as a sum of N terms different terms, 
H = Y^=i rms H i, where ||#;|| < 1 for all i. Given any density matrix p, we define the 
energy density of p to be 

trjpH) - E 

N terms ' (8) 

where E is the lowest eigenvalue of H. 

Note that the energy density may depend upon the decomposition used to write H. 

A simple way to construct ipi ow is as follows. Take the L-hj-L square lattice and divide 
it into small squares of size /-by-/, for / << L (if L/l is not an integer, we can use a 
mix of squares and rectangles of size / and / + 1). Label the small squares by an index a 
and consider the state p = Y[ a P<n where p a is a mixed state: it is the equal probability 
mixture of the state with all spins up in square a and the state with all spins down, so 
p a = (1/2) (I ttt — )(ttt —I + I III •••)(!!! •••!)• Then, this state p has an energy density 
of order l/L which can be made arbitrarily small for large L, but the correlation function 
tr(pS?Sj) is equal to is sites i, j are in different small squares. 

We could define long-range order in this case in a slightly different way. Consider the 
quantity 

m2 = ^Gl^ 2 ) 2 ) 2 ' ( 9 ) 

i 

where N = L 2 is the total number of spins. Using the checkerboard state defined above, we 
have M 2 = if L/l is an even integer. For the mixed state n<tPa> we nave ^ 2 ~ I 2 / L 2 « 1. 
Note that this quantity M 2 is equal to the average, over all choices of i,j, of the correlation 
function tr^pS^S*). On expanders (below), M behaves quite differently. 

Note that nothing in the above discussion of the checkerboard state or state Y[ a Pa relied 
crucially on having a two-dimensional lattice. We could have made a similar construction 
using small hypercubes of linear size / for any finite dimensional lattice. 

So, finally, we consider an Ising model on an expander graph. We define the expansion 
coefficient by 

Definition 3.2. Suppose the graph has N vertices and has the expansion property that for 
every set S of vertices with \S\ < N/2, the number of edges in the graph connecting S to 
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its complement has cardinality at least c\S\ for some constant c > 0. Then the graph has 
expansion coefficient c. 

Suppose c > 0. Then, we claim for any state of sufficiently low energy, the quantity M 2 
is of order unity, in contrast the the case on the square lattice. Consider a configuration 
of spins with spins up (that is, with Sf = +1) and N± spins down with + N± = N. 
Assume that iVj- > N±. Then the set of vertices with down spins has cardinality at most 
\S\/2 and hence the number of edges connecting the down spins to the up spins is at least 
a/Vj_. In general, we find that the number of edges connecting down spins to up spins is at 
least min(A^, N±). Let the energy of the ground states of the system (either the state with 
all spins up or the state with all spins down) be E . Note that E is equal to —1 times the 
number of edges. Every edge connecting an up spin to a down spin raises the energy by 2 
compared to the ground state energy so the energy E of a configurations with up spins 
and iVj, down spins obeys 

E > E + 2cmin(iV t , N±). (10) 



Note that 



M*=W^£. (11) 
Assuming without loss of generality that iVj. < A^, we find that 

M 2 = (N-2N l ) 2 /N 2 (12) 

> (N 2 — 4NN±) /N 2 

= 1 - 4NJN 

2(E - E ) 

cN 

For a mixed state p, we can average both sides of the above equation and find that 

M*>l-?^=^, (13) 

where now E = tr(pH Ferro ). Note that the number of terms in H Ferro is equal to (N/2) 
times the degree d of the graph. So, (E — E )/N is equal to the energy density times 2/d, 
and so indeed states with small enough energy density have M 2 of order unity. 

We now turn to the case of topological order, considering the toric code on various 
lattices. First consider the toric code on a two dimensional lattice. Depending on boundary 



21 



conditions chosen (periodic, i.e. torus boundary conditions, open boundaries, etc..) there 
may be one or more zero energy ground states. However, these states are nontrivial under 
the circuit definition, in the sense that none of them is (R, e)-trivial for R/L sufficiently 
small and for e sufficiently small compared to unity. See for example Ref. [2 for a torus 
geometry or Ref. |9| for other geometries. However, at any fixed nonzero temperature the 
states are trivial. More precisely, for any such temperature, for any L, the state is (R, 1/2) 
trivial for an R which grows only logarithmically with L (the prefactor in this logarithmic 
dependence depends upon temperature). See Ref. lo|, for example. 

Turning to the case of the four- dimensional lattice, there are various possible toric codes 
that can be defined. In general, in d dimensions, we can refer to a (q, r) code, with q + r = d, 
where q — 1 is the dimensionality of electric defects and r — 1 is the dimensionality of magnetic 
defects. The electric defects bound g-dimensional surfaces, and the magnetic defects bound 
r-dimensional surfaces. The toric code discussed above in two dimensions is an example of 
a (1, 1) code. In the next section, we describe a generalization of this toric code to arbitrary 
two-complexes. However, for now we consider a (2,2) code on the four-dimensional lattice. 
In this case, for sufficiently small nonzero temperature, the state is not (R, e) trivial for R/L 
sufficiently small and e sufficiently small compared to unity^ 1 ^. 

So, we have seen three classical systems, one of which has long-range order only at zero 
temperature, the next has long-range order at sufficiently small temperature, and the third 
has order for all sufficiently low energy states. We have quantum systems that are in a sense 
analogous to the first two of these. The question naturally arises then: is there a quantum 
system that has NLTS? This is the question we address. 

4. TORIC CODES ON COMPLEXES 

For any 2-complex, we can define a Hamiltonian that generalizes the toric code Hamilto- 
nian as follows. Associated with every 1-cell of the complex is a spin-1/2 degree of freedom. 
The Hamiltonian is a sum of terms 

H tc = -A-B, (14) 

where A, B are defined by 



A = -Y,A S , (15) 
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B = -J2B P . (16) 

v 

The sum in Eq. (fl5]) is a sum over all 0-cells, where s labels the 0-cell, and A s is defined to 
by 

A s = n ( i7 ) 

sGde 

with the product being a product over all 1-cells e that are attached to s (so that s is in the 
boundary de of e) and 5^ is the spin operator S x acting on the degree of freedom on e. The 
sum in Eq. f fl6|) is a sum over all 2-cells, where p labels the 2-cell, and B p is defined to by 

b p = n s p> 

with the product being a product over all 1-cells e that are attached to p and S% is the spin 
operator S z acting on the degree of freedom on e. We have [A s , B p ] = for all s,p, for all 
2-complexes. 

We refer to the terms A s as vertex terms and the terms B p as plaquette terms. If there 
is a state ip and a 0-cell s for which A s ip = —ip, we say that if) has a vertex defect. The 
number of vertex defects of ip is given by Ylsi^i — A s )i/))/2. We use similar terminology 
for plaquette defects. 

The construction above is given for an arbitrary complex. If the complex is a triangulation 
of a manifold, then the above code is a (1, d — 1) code. 

5. TORIC CODE ON G x G: TOWARD THE QUANTUM PCP CONJECTURE 

We now consider this Hamiltonian (TH1) on the 2-complex obtained from G x G defined 
previously. The graphs G will all have bounded degree and diverging girth. We will prove 
a "one-sided" result on the absence of low energy trivial states with no vertex defects, 
theorem 15.41 Later, in section l5Al we discuss a possible alternative proof idea of our result; 
this alternative proof idea relies on an unproven conjecture but may be interesting as, if this 
conjecture is proven, this method may be useful for a variety of other Hamiltonians. 

Definition 5.1. Give a set S of Hermitian operators which square to the identity, we say 
that a state ip is stabilized by S if Pijj = ip for all operators P in S. 
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Note that in this definition that S may be an arbitrary set of such operators. There 
is a formalism called the stabilizer formalism based on operators which are products of 
Pauli operators, but we do not make this restriction here. We also do not require that the 
operators commute with each other. 

Suppose ip is a trivial state. Then, ip = U C i rcu iti for some product state £ and for some 
unitary U C i rcu u defined by a local quantum circuit. The product state £ is stabilized by a set 
of operators Oi, one such projector for every degree of freedom i. Each such operator Oi is a 
tensor product of an operator on degree of freedom i with the identity operator elsewhere; the 
operator on site i has a unique eigenvalue equal to +1, so that specifying that ip is stabilized 
by the set of operators Oj uniquely specifies ip U P to a scalar. In the case of the toric code, 
since the degrees of freedom are on 1-cells, the index i ranges over all 1-cells. Therefore, 
ip is stabilized by the set of operators V T C ircuit0^j\ ircuit . Define Oi = U C i rcu i t OiUl ircuit . Note 
that if the quantum circuit has range R, then each operator O, has support on the degrees 
of freedom within distance R of i. If ip minimizes A, then A s ip = ip for all s, and so ip is 
stabilized also by the set of operators A s . 

We now use this idea of stabilizers to define a new state <p in lemma 15.31 Very roughly 
speaking, <p will have a strange property: in local regions it will "look like" ip does (that 
is, have the same reduced density matrix), up to possible conjugation by a certain unitary 
operator that anti-commutes with certain plaquette operators B p . Thus, the unitary op- 
erator can be viewed as "inserting defects" on the plaquettes that it anti-commutes with. 
However, there will be a seemingly contradictory property. To explain this property, note 
that we are used to the requirement that defects be inserted in pairs in the two-dimensional 
toric code on a torus, and every state must have an even number of plaquette defects (this 
is a consequence of the second homology of the torus), and in this G x G there is a similar 
requirement: due to the branching structure of G x G, things are more complicated but still 
not all configurations of defects are possible. However, the conjugation by unitary that we 
use to define <p will appear not to respect this property: it will seem to add defects in certain 
places but not add corresponding defects in other places. The resolution of this apparent 
paradox will turn out to be that the original state ip had defects, and in fact we will use 
this construction of (p to prove that ip has defects, and indeed to prove a lower bound on 
the energy density of ip for the toric code on G x G. 

Before defining the state 0, let us define the needed unitary operators. 
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Definition 5.2. Let p be any 1-cell in G x G. Define the operator D x {p) as follows (this 
construction involves an arbitrary choice of a 1-cell below; this choice can be made in any 
way, but every time we refer to the operator D x {p) for given p, we take the same choice). 
A 1-cell in G x G is labeled by an edge in one graph and a vertex in another graph. Pick 
a 1-cell attached to p which is labeled by an edge a in the first graph and a vertex i in the 
second graph. So, we refer to this 1-cell as (a,i). Pick another 1-cell attached to p which is 
labeled by an edge b in the second graph and a vertex j in the first graph. So, we refer to 
this 1-cell as (j,b). Note that edge a is attached to vertex j in G and edge b is attached to 
vertex i in G. 

If the graph G has degree d, then i has d — 1 edges attached to it other than b. Let R be 
the set of all vertices in G which are within distance g/2 of i, where g is the girth of G, such 
that the shortest path from i to the given vertex does not use edge b. There are (D — I) 9 / 2 
such vertices in R. 

Define the operator D x {p) to be the product of a x over all 1-cells e of the form (a, k) for 
k in R. 

We define another operator C x {p) as follows. Define T to be the set of all vertices in G 
which are within distance g/2 of j , such that the shortest path from j to the given vertex 
does not use edge a. Let OR be the boundary of R: the set of edges in G which join vertices 
in R to vertices not in R; let dT be similarly defined as the boundary of T . Define C x (p) 
to be the product of a x over all 1-cells e of the form (k,b) for k in T and over all 1-cells e 
of the form (c, k) for k in R and c in dT and over all 1-cells e of the form {k, c) for k inT 
and c in dR. 

We note that the product C x {p)D x {p) is equal to a product of operators A s over an 
appropriately chosen set of sites s. That, the two 1-chains corresponding to these two 
different operators are homologous. Further, C x {p) anti-commutes with B p , as well as with 
certain other plaquette operators a distance g/2 from p (the plaquette operators it anti- 
commutes with are those in the coboundary of the 1-chain corresponding to C x (p)). It 
commutes with all other B p . The operator D x (p) has the same commutation and anti- 
commutation relations with B p . 

These two different operators C x (p), D x (p) can be thought of in the following way: con- 
sider the two dimensional toric code. A product of a x along a line of bonds on the dual 
lattice creates a pair of anyons at the ends of the line (i.e., thinking of the line as a 1-chain, 
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it creates anyons in the coboundary); there are various possible choices of lines homologous 
to the given line, and for each such line we can define an operator which is the product of o~ x 
along that line. The operators C x (p) and D x (p) are like two such choices of lines, however 
due to the branching nature of the lattice, rather than creating a pair of defects they have 
many plaquettes in their coboundary. 

Lemma 5.3. Let if) be a trivial state stabilized by the set of operators Oi as defined above. 
Let ip also be stabilized by the set of operators A s . Let V be any set of plaquettes. Assume 
that the range R of the quantum circuit is less than g/32. Then, there exists a state <f>, with 
\<f>\ = 1, such that <f> is stabilized by the set of operators A s and also is stabilized by the set 
of operators 

( n Dx (p))Oi( n ( i9 ) 

pe"P,dist(p.i)<#/4 pe7 , ,dist(p,i)<g/4 

where the notation on the product means the product over allp in V such that dist(p, i) < g/4. 
Note that the operator (j^[ peV dist ( p i)< g /4 D x {p) s j is a Hermitian unitary. 

Further, consider any set X of 1- cells inG x G, with diameter of X at most g/8. Let Vx 
be the set of plaquettes in V such that all 1-cells attached to the plaquette are in X. Let (f)x 
be the reduced density matrix of the state <f> on X and let if>x be the reduced density matrix 
of the state if) on X . Then, 

= ( n dx (p))^( n d *(p))- (2°) 

P&Vx P&Vx 

Proof. We show Eq. f[2TJ|) . showing the equation under the assumption that <f> exists. Then 
we prove the existence of <f>. 

First, we note an important property. There is an arbitrariness in how we chose the 
stabilizers. For each i, define an arbitrary set of 0-cells called Si. Then, if we specified that 
a state is stabilized by the set of operators A s and also is stabilized by the set of operators 

(n^)( n Dx (p))Oi( n ^>))(n^)' (2i) 

seSi p£V,dist(p,i)<g/4 pe7 , ,dist(p,i)<g/4 sESi 

then it is also specified by the operators in Eq. (fl9l) and conversely if it is specified by 
the set of operators A s and also is stabilized by the set of operators in Eq. fTTTJT) then it 
is also stabilized by the set of operators in Eq. fT2TT) . This can be regarded as a "gauge 
freedom". Since, as we have noted before, the product C x (p)D x (p) is a product of A s over 
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an appropriate set of s, this means that we can replace operators D x {p) in Eq. ([19]) by 
operators C x (p) as desired, so long as we replace all occurences of D x (p) for a given i,p. 
That is, for each z", we can define a set T>i with T>i C V and then the state can be taken to 
be stabilized by the operators 



( n n cx (p))Oi( n ° x ^){ n ° x ^)- 

pe£>i,dist(p,i)<g/4 peV\Vi,dist(p,i)<g/4 peV\D i ,dist(p,i)<g/'i p£Vi,dist(p,i)<g/4 

(22) 



Consider a given z. Note that if dist(p, i) < g/A and the operator Oi is supported on a 
set of diameter at most R, then one of the following cases holds. Either, [D x (p),Oi\ = or 
z is within distance R of the support of D x (p). However, this latter case means that either 
% is not within distance R of the support of C x (p) or i is within distance R of p. Thus, 
there are three possible cases if dist(p, z) < g/A: Oj commutes with D x (p), or Oj commutes 
with C x (p), or i is within distance R of p. If Oi commutes with D x (p) then we can remove 
all occurences of D x (p) from Eq. (1221) . If Oj commutes with C x (p), then we can conjugate 
Eq. (1221) by C x (p), and then conjugate by a product of operators A s to turn that C x {p) into 
a D x {p). Thus, for any z, we can define any set of plaquettes V{ such that every plaquette p 
in V% has dist(p, z) < g/A and such that every plaquette pmV with dist(p, z) < i? is in Pj, 
and then the state is stabilized by the set of v4 s and by the set of operators 



This is the final form of the stabilizers that we will take. One way to describe these various 
choices of stabilizers is to note that the stabilizers generate a group and every element of the 
group stabilizes ip: these different choices of the stabilizers correspond to conjugating certain 
generators of the group by other generators of the group, giving a different presentation of 
the same group. However, now Eq. (1201) follows almost immediately. For each i such that 
dist(z", X) < R, take V\ to include all plaquettes p G V within distance 2R of X. Thus, 
Vi has the needed property of including all plaquettes within distance R of i. Also, by a 
triangle inequality, since z is within distance R of X, every plaquette in Vi is within distance 
3R + g/8 of z, and for the given choice of R this is less than g/A. 

We now prove that such a state <fi exists. The proof is inductive. The induction is on the 
cardinality of V . It is clear that the claim holds when V is the empty set, since then we can 
take 4> — ij}. So, we must prove that if the claim holds for a given set V, then it also holds 




(23) 
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for the set of plaquettes V = V U {q}, for any plaquette q. We now show this. Assume that 
such a <p exists for a given set V. Thus, there is a state stabilized by the set of operators A s 
and by the set of operators 

( n Dx (p))Oi( n ^(p))- ( 24 ) 

pG'P,dist(p,i)<g/4 peP,dist(p,i)<g/4 

Hence there is a state stabilized by the operators A s and by the set of operators in Eq. f )24|) 
for dist(g, i) > R. (This last statement is a triviality: since it is stabilized by the operators 
of Eq. ( 124"|) . it is stabilized by a subset of those operators). 

Note, then, that there is a state stabilized by the set of operators A s and by the set of 
operators 

( n Dx (p))Oi( n d ») ( 2s ) 

p€'P / ,dist(p,i)<g/4 p£~P' ,dist(p,i)<g/4 

for dist(g, i) > R. Here we have used the fact that for dist(g, i) > R and dist(g, i) < g/4, we 
can choose to conjugate a given stabilizer by the operator D x (p) (or choose not to conjugate 
it) without changing the group of stabilizers. Similar to before, we can define a set of 
plaquettes V[ for each i, such that every plaquette p in V[ has dist(p, %) < g/A and such that 
every plaquette p in V with dist(p, i) < R is in V[ such that the state is stabilized by the 
set of A s and by the set of operators 

(n^(p))^(n^)) ( 26 ) 

for dist(g, i) > R. For each i such that dist(g, i) < g/8, take V[ to be the set of all plaquettes 
p EV' within distance R + g/8 of q. 

Now take this set of stabilizers Eq. ( |2"oT) and conjugate all stabilizers by ripgp' D x (p) for j 
chosen to be any 1-cell with dist(g, j) < g/8 (note that all such 1-cells give the same set Pj). 
This new set of stabilizers will stabilize a different state from the state stabilized by Eq. (|26|) : 
it will stabilize that state multiplied by ripeT" D x {p). The resulting state is stabilized by 
the set of A s and by Oi for R < dist(g, i) < g/8 and also by various other stabilizers. These 
other stabilizers are obtained by conjugating the stabilizers for dist(g, i) > g/8. We do not 
bother writing them down, but note the important fact that the support of those stabilizers 
for any given i is the same as the support of Oi. Take all these stabilizers and conjugate 
them by U\ ircuit ] since U circui t need not be Hermitian, let us specify that by "conjugate" we 
mean left multiply by U\ ircuit and right multiply by U circuit- The resulting state is stabilized 
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by operators U circuit A s U C i rcuitl by operators Oj for R < dist(g, i) < g/8, and by various other 
operators whose support is a distance at least g/8 — 2R away from q. Now, if we could show 
that a state existed that was stabilized by the same operators and also by the operators Oi 
for dist (g, i) < R then we would be done: we simply undo the conjugation of the stabilizers 
by U\ ircuit and by ripcP' D x (p) and we have a state that is stabilizers by the desired set 
of stabilizers. However, we claim that such a state does exist: first, we begin by removing 
some of the stabilizers: we require that the state be stabilized by Ul ircuit A s U circui t only for 
dist(s,g) > g/16, while continuing to require that it also be stabilized by operators Oi for 
R < dist(g, i) < g/8, and by various other operators whose support is a distance at least 
g/8 — 2R away from q. Note now that all the operators Ul ircuit A s U circuit have support at 
least a distance g/16 — R > g/32 away from q. Therefore, the only stabilizers which have 
support within g/32 of q are the operators Oi. Each of these operators acts on a single 
1-cell. Thus, no stabilizer has support on the 1-cells i for dist(i,g) < R. Thus, we can 
also require that the state be stabilized by the operators Oi for dist(g, i) < R. Thus, we 
have a state stabilized by Ul ircuit A s U C i rcu i t only for dist(s,g) > g/16, and stabilized by Oi 
for dist(g, i) < g/8, and by various other operators whose support is a distance at least 
g/8 — 2R away from q. Undo the conjugation by Ul ircuit and by n p ep' D x {p). This gives a 
state stabilized by A s for dist(s,g) > g/16 and also conjugated by the set of operators in 
Eq. (125]) for all i. We simply need to show that the state is also stabilized by the A s for 
dist(s, q) < g/16. However, note that being stabilized by the set of operators in Eq. ( |25i) for 
i with dist(z, q) < g/16 + R + 1 uniquely specifies the reduced density matrix of the state 
on the set of 1-cells within distance g/16 + 1 of q and hence on the set of 1-cells attached 
to all 0-cells within distance g/16 of q. Further, we know that a state does exist which is 
stabilized by the set of operators in Eq. (1251) for i with dist(i, q) < g/16 + R + 1, as such a 
state can be obtained by unitarily conjugating %p, and this state is stabilized by A s for all s 
within distance g/16 of q. Hence, the state that we have constructed is indeed stabilized by 
all the A s . □ 

We now show that there are no low energy trivial states with no vertex defects. 

Theorem 5.4. There is an e > such that for all R < oo, for all sufficiently large graphs 
G in the family of graphs defined above there is no R-trivial state ip on G x G such that 
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A s ip = ip for all ip and = 1 and such that 

(^-B^)<-(l-e)N p , (27) 
where N p is the number of plaquettes. 

Proof. Pick a set V of plaquettes which will soon be treated as a variable. For given R, for 
sufficiently large G, define the state as in lemma 15.31 Define B' by 

B> = -Y,B p + J2 B p- ( 28 ) 

That is, B' is also a sum of plaquette terms, just like B is, except the sign of the terms is 
flipped for plaquettes in V . Note that by Eq. (120]) . 

(V,-SV») = ^,S / 0>. (29) 

However, we claim that since b2{G x G) is extensive, there is an e > such that for all 
sufficiently large G we can find a set of plaquettes V such that the minimum eigenvalue of 
B' is greater than or equal to — (1 — e)N p . Once we have proven this claim, this will complete 
the proof, since it will imply, by Eq. ( 1291) that 

(i/>,-B1>)>-(l-e)N p (30) 

for all sufficiently large G. 

To prove the claim on the minimum eigenvalue of B', note that B' is diagonal in the 
S z basis so every eigenvector of B' can be chosen to be an eigenvector of all the operators 
SI simultaneously. Let N e be the number of 1-cells. Then, there are 2 6 different states 
which are eigenvectors of all the S* separately: these states correspond to the 2 Ne different 
assignments of +1 or —1 to each 1-cell. These states form an orthonormal basis. 

Each such basis state £ is an eigenstate of all the B p . Number the plaquettes p from 
l...N p . So, to each such basis state £ we can define a vector with N p different entries, 
with the p-th entry of the vector being +1 or —1 depending upon whether the expectation 
value of B p is +1 or —1. Each choice of a set V can be used to define another vector v-p: 
we define this vector to be +1 on the p-th plaquette if p G V and —1 otherwise. Then, the 
expectation value 

(£,B'0 = -N p + J2h-Vv\. (31) 

p 
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The term on the right-hand side of the above equation can be understood also in terms of a 
Hamming distance. Given a vector with entries +1 or —1 we can replace each +1 with a 
and each —1 with a 1 and the result is a bit string. Then, the term on the right-hand side of 
the above equation is equal to twice the Hamming distance between the corresponding bit 
string. For a given basis state £, let be the corresponding bit string; i.e., the bit string 
obtained from v% by replacing +1 with and —1 with 1 and let sp be the bit string obtained 
from v-p by a similar replacement. 

There are 2 Np possible vectors v-p. On the other hand, there are at most 2 Ne possible 
vectors v$ (in fact, the number is slightly smaller than this since the mapping from basis 
states £ to vectors v% is not one-to-one). So, since N p > N e , there certainly are choices of V 
such that the distance on the right-hand side of the above equation is positive; i.e., choices 
of V such that the minimum eigenvalue of B' is greater than —N p . However, we need a 
stronger result than this. This stronger result follows by a similar counting argument. For 
every bit string of length N p , the number of bit strings within Hamming distance (e/2)N p 
of the given bit string is roughly (up to subexponential corrections) 

exp(N p H(e/2)), (32) 

were H(x) is the entropy function: H(x) = —xln(x) — (1 — x) ln(l — x). Thus, the number 
of bit strings within Hamming distance N p of the set of string is at most (up to subexpo- 
nential corrections) 2 Ne exp(N p H(e/2)). For sufficiently small e, since N p /N e is greater than 
one, we can ensure that 2 Ne exp(N e H (e/2)) < 2 Np . Thus, for sufficiently small e, there is 
some V such that s-p has Hamming distance greater than e/2 from all s$. Thus, for this V, 
for all basis states £ the expectation value (£, B'£) is greater than —(1 — e)N p . □ 



A. An Idea for An Alternate Proof Using Expansion Properties 

Let S be the coboundary operator. Suppose the following conjecture holds: "efficient 
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Conjecture 5.5. There is a constant c > such that the following holds for all sufficiently 
large graphs G form the family considered above. Let x be any 1-chain with Z 2 coefficients 
on G x G. Then, there is a 1-chain y which is a cocycle such that that 

|supp(<5x)| 



|supp(a; + y)| 
31 



> c, (33) 



where 5 is the coboundary operator, supp(...) is the support of the coboundary of a vector 
and |supp(...)| denotes the cardinality of the support of a vector. 
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but it is 



This conjecture is similar to a cohomological expansion definition in Ref. 
crucially different in that we allow y to be any cocycle rather than just a coboundary 
Suppose this conjecture holds. Under this assumption, we sketch an alternative proof of 
theorem 15.41 Since the conjecture itself is unproven, we will given only a sketch and many 
details will be left out. 

To motivate the alternate argument, suppose first that we had a state tp such that B p ip = 
ijj for all p. Then, there would be a certain kind of long-range order in the state ip. Consider 
a "Wilson loop" operator, W(C), which we define to be the product of S* around a 1-cycle 
C. If the 1-cycle is contractible, then (if), W(C)ip) = +1. Consider instead cycles C which 
are the product of a given vertex v in the first graph with a cycle c in the second graph. 
We will write C(v, c) to denote such a cycle, determined by the given vertex v and cycle c. 
Then, 

y>,W(v,c)4>) = (,J>,W(v',cW) (34) 

for any pair of vertices v,v' and any cycle c. The first step of our sketch is to show that 
any state with low energy the plaquette terms has a similar (but weaker) kind of long-range 
order, defined below. 

Let x be a 1-chain. Consider first a state ip c i{x) which is a product state, which is an 
eigenstate of all operators with S% having expectation value —1 for e in the chain and 
+1 otherwise. Here the subscript cl stands for "classical"; we refer later to these states 
as classical states and they form an orthonormal basis. Suppose that (ip c i(x), —Bip c i(x)) < 
— (1 — e)N p . Then, by the conjecture above, we can find a y that is a cocycle such that x + y 
has a small support. That is |supp(x + y)\ < eN p /c. Consider the state tp c i(y)- This state 
obeys B p ip c i(y) = ip c i{y) for all p, so it has long-range order in the Wilson loop operators as 
in Eq. AM]) . 

We can restate the long-range order in Wilson loops in a slightly different way. Let the 
graph G have N edges and degree d. Hence the graph has first Betti number b\{G) = ^=2 + 1. 
For each vertex v in G, we define a vector s(v) with _|_ \ entries; these entries are the 
expectation values (ip c i(x), W(v, c)ip c i(x)) for a set of cycles c that form a basis for the first 
homology. The entries in this vector equal ±1. 
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Note that if we instead consider a vector t(v) with entries being the expectation values 
)ip c i(y),W(vic)ipci(y)), then t(v) = t(v') for all v,v'. The chain x + y has small support. Let 
N(x + y, v) denote the number of 1-cells in x + y which are a product of the given vertex v 
in the first graph and an arbitrary edge in the second graph. Note that Yl e + Vi v ) — 
|supp(x + y)\. So, speaking loosely, for most v, we have N(x + y,v) small compared to N. 
For example, for at least 1/2 of v, we have 

N( X + y )<2^ = 2 C -^N. (35) 

Note that N p /N 2 = 0(1), so for at least 1/2 of v, N(x + y,v) is bounded by 0(l)eN = 
0(l)eN e . (The specific number 1/2 is not important.) 

We now define a distance function and use this bound on N(x+y, v) to bound the distance 
between s(v) and s(v'). Note that the vector s(v) is completely determined by the vector of 
expectation values (ipd, S*?p c i) for 1-cells e which are a product of the given vertex v in the 
first graph and an arbitrary edge in the second graph. This second vector has N e entries. 
In fact, there are many vectors of expectation values (i/j c i, S*ijj c i) which determine the same 
vector s(v). In general, we will say that u s(v) is the Wilson loop vector of w" if the entry 
of s(v) corresponding to cycle c is equal to the product of the entries of w corresponding to 
the edges in that cycle. We define a distance function between two vectors, s(v) and s(v') as 
follows. Take the minimum, over all vectors w,w' such that s(v) is the Wilson loop vector 
of w and s(v') is the Wilson loop vector of w', of the Hamming distance between w and w' . 

Then, we have the for at least 1/2 of v, the distance between s(v) and t(v) is at most 
0(l)eN e . Then, for at least 1/4 of the pairs v,v' the distance between s(v) and s(v') is at 
most 20(l)eN e = 0(l)eN e . The number of possible choices of at most 0(l)eiV e edges from 
a set of iV e edges is exp(0(l)eA^ e ), and for at least 1/4 of the pairs v, v\ given s(v) the vector 
s(v') is one of at most exp(0(l)eA^ e ) possible vectors. 

We now define a quantum operator that we write as 8w(v),z)- The notation is intended to 
be suggestive here: z is an arbitrary vector with bi(G) entries, the entries equaling ±1, and 
this operator will compute the probability that the vector of Wilson loop expectation value 
is equal to z. We define 

tt l + W{v,c)z c 
o W (v),z = H g t 36 ) 

c 

where the product is over a set of cycles c that form a basis for the first homology and z c 
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denotes the entry of z correspnding to cycle c. The operators 

1 + W(v,c)z c 
2 

are projectors which commute with each other, so 5w(v),z is a projector. Note that 



(37) 



(ipd(x),S w{v)tZ i/j c i(x)) (38) 

equal 1 if z — s(v) and otherwise. 

Compute the sum of expectation values 

dist(z,z')<0(l)eN e 

Y {Mx),6 W ( iV ) ig 6 a ( v >) >g >il>d(x)) (39) 

z,z' 

where the sum is over vectors z, z' with bi(G) entries each having values ±1. For randomly 
chosen v,v', this expectation value is at least 1/4. This is the "long-range order" that we 
identify for this state ip c i. We now consider any arbitrary state ip (not necessarily a classical 
state) such that < —(1 — e/2)N p . This state ip is a superposition of classical 

states: ip = J2 X A(x)ip c i(x) . Choosing a classical state ip c i( x ) with probability distribution 
|A(x)| 2 , with probability at least 1/2 we find that (ip c i(x), — Bip d (x)) < —(1 — e)N p Thus, 
for any such state ip, we find that the expecation value 

dist(z,z')<0(l)eN e 

Y (^,Sw(v),zS s (v'),z'^) (40) 

z,z' 

is at least 1/8. 

Note, however, that if tp is an i?-trivial state, then the expectation value factors if the 
distance between v and v' is greater than 2R (note that here by distance we are referring 
to the distance on the graph, using the graph metric). Further, for a randomly chosen pair 
v, v', the probability that the distance between v, v' is greater than 2R tends to 1 as iV — > oo. 
Thus, up to corrections that go to zero as iV — > oo, we find that 

dist(z,z')<0(l)eN e dist(z,z')<0(l)eN e 

S W(v),zS s (v'),z'^) = Yl ( 41 ) 

z,z' z,z' 

Finally, suppose that ip had the property that for all v we had that 

(VS<W^ = 2 ^ 1(G) - ( 42 ) 
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That is, suppose that there was a uniform probability distribution, independent of z, of 
having different Wilson loop vectors. This can be seen to contradict Eq. (j41j) for sufficiently 
small e, since in that case the fraction 

exp(Q(l)eiV e ) 

26i(G) 1 J 

tends to zero as iV e — > 00 for small enough e and this fraction represents the fraction 
of vectors within distance 0(l)eN e of a given vector. Hence, we cannot have a uniform 
probability distribution. 

However, it can be shown that having a non-uniform probability distribution contradicts 
the assumption of a trivial state with no vertex defects. The argument involves constructing 
operators similar to the C x (p) and D x (p) above. For any 1-cell e in G x G which is a product 
of a vertex v in the first graph and an edge a in the second graph, we can construct an 
operator O that anti-commutes with S* and which commutes with all other Si for / being a 
product of vertex v with an edge b 7^ a. If the girth of the graph is large enough compared to 
R, we can construct this operator O such that it commutes with all the plaquette operators 
Bp except for certain plaquttes p far from vxG (note: this means that the plaquettes are far 
from all 1-cells e of the form v times an edge in the graph, not just far from e); this operator 
O will anti-commute with B p for p in some set S of plaquettes with dist(S', v x G) > 2R 
and will commute with all other B p . Just like the C x (p) and D x (p) were, this operator O 
will be a product of SJ over / in some set of 1-cells. We can also construct an operator 
O' such that 00' is a product of operators A s . Then, 00' ip = ip. This argument that we 



Further, we can choose O, O' 



now follow is quite similar to that in section VIII of Ref. 
such that the intersection of the support of O with the set of sites within distance 2R of 
v x G is at least distance 2R away from the intersection of the support of O' with the set 
of sites within distance 2R of v x G. Using the assumption that the state is i?-trivial, we 
can show that ip and Otp have the same reduced density matrices onoxG (to see this, 
pull back on the quantum circuit, so that we consider operators O = Ul ircuit OU C i rcuit and 
O' = Ul ircuit O'U circuit acting on a product state ; using the distance between the supports of 
O and O' mentioned above, we can see that on the set of 1-cells in v x G), the supports of O 
do not overlap; since 00' acting on the product state gives the product state, this implies 
that O leaves the reduced density matrix on«xG unchanged). So, 

(f.^M^) = (ip,OS w ^ v)>z )Oip) = (ip,5 W (v),z4), (44) 
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where z' is the vector obtained by finding a w such that z is the Wilson loop vector of w, 
then flipping the value of the entry in w corresponding to edge a to obtain a vector w' and 
then letting z' be the Wilson loop vector of w'. Using the above equation, for various choices 
of a, we can see that (ip, SwM^ip) indeed must be independent of z for this .R-trivial state. 



B. Triviality of Low Energy States with Vertex Defects 

Theorem 15.41 showed the absence of trivial low energy states without any vertex defects. 
Now, we show the converse: there are trivial low energy states with vertex defects. Indeed, 
there are trivial states with any arbitrarily low density of vertex defects and with no plaquette 
defects. 

Theorem 5.6. For any e > 0, there exists an R such that for all large sufficiently graphs 
G in the family of graphs defined above there is an R-trivial state ip on G x G such that 
B p ip = tp for all ij) and such that 

(^-A^)<-{l-e)N a> (45) 
where N s is the number of plaquettes. 

Proof. We construct the state ip by picking a subset S of the vertices, with the subset 
containing a fraction e of the vertices, and construct an .R-trivial state that is an exact 
ground state of the Hamiltonian 

H tc {S) = -B~Y J A s . (46) 

That is, it is an exact ground state of the Hamiltonian with the operators A s for s G 5 
removed. 

We construct this R- trivial state by giving a unitary quantum circuit that turns H tc (S) 
into a Hamtilonian that is diagonal in the S z basis and hence which clearly has a product 
state ground state. The procedure that follows is closely related to that in Appendix A of 



Ref. |10j. Consider a given vertex s £ S. Consider any 1-cell e attached to this vertex. Note 
that there are two 0-cells attached to e and so there is some other 0-cell s' attached to e 
with s' 7^ s. Assume that s' (jL S. Define the unitary transformation 

U e , s , = e W {\s z e A s ,). (47) 
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Then, 

U e , s Htc(S)Ul s , = H tc (S U {s'}) - S z e . (48) 

Define a unitary U s to be the product of C/ eiS /( e , s ) over a subset (given in the next sentence) 
of the 1-cells e attached to s with s'(e, s) being the 0-cell other than s that is attached to e. 
The particular subset we choose is the subset of e such that s'(e, s) S. 

U s H tc (S)Ul = H tc (S'), (49) 

where S' is equal to the union of S with the set of all 0-cells attached to a 1-cell attached 
to s. 

We now define the circuit. Let S = S. Let H = H(S). Let Hi = UiHi-iUj. We will 
choose the unitaries Ui so that each Hamiltonain Hi is equal to H tc (Si) plus some sum of 
operators SI, with the sum being taken over a set of edges Ei. The set of edges Ei will 
have the property that every edge in Ei is not in the support of any operator in the sum 
— Y2sgsAs, and hence the Hamiltonian will still be a set of commuting terms. 

For given e, we pick S to be a set of s such that every s' is within some bounded (e- 
dependent) distance of So- The i-th round of the circuit will be a product of U s over some 
given set Xi of s. We choose a set X^ such that there is no pair s, s' G X with dist(s, s') < 2. 
We claim that we can choose such sets Xi such that after a bounded number of rounds, all 
0-cells are in S, leaving a Hamiltonian that is diagonal in the S z basis. Here is one way 
to show this claim: color the set of all 0-cells by some bounded number of different colors 
such that no two 0-cells within distance two of each other have the same color. This can be 
done with c = 0(1) colors. Let the i-th round of the circuit have Xi all i with color i mode 
such that s G 5,. Then, after c rounds, S c will contain all s within distance 1 of So, after 
2c rounds, S2 C will contain all s within distance 2 of So, and so on, and so after a bounded 
number k of rounds, Sk will contain all s. □ 



6. HYPERGRAPH CODES 



This above section naturally leads to the question of whether we can construct a Hamil- 
tonian which has NLTS, without making this "one-sided" assumption. One interesting way 
to attack this question is to consider the hypergraph product codes of Ref. [l]. We will use 
a slightly different language, the language of chain complexes, to define these codes. This 
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language will make certain generalizations possible, as we discuss in the next two subsec- 
tions. Ref. [l uses the language of hypergraphs; we begin by describing hypergraphs as chain 
complexes. 



A. Hypergraphs and Hypercomplexes 

We can define a hypergraph with hypervertices and hyperedges by a bipartite graph. The 
graph is divided into left and right parts, with the left part of the graph representing the 
hyperedges of the hypergraph and the right part of the graph representing the hyperver- 
tices, with an edge in the graph between a right node and a left node if the hypervertex 
corresponding to the right node is in the hyperedge corresponding to the left node. Given a 
graph representation of a hypergraph, we can define a space C\ (this space may be a vector 
space or it may be more general; see below) with dimensionality equal to the number of left 
nodes and a space Cq with dimensionality equal to the number of right nodes, and a linear 
operator d\ called the "boundary operator" from G\ to C , with the linear operator di being 
the adjacency matrix of the graph. This representation suggests a natural generalization of 
hypergraphs to "hypercomplexes" . Such a hypercomplex in fact is simply a chain complex 
with certain sparseness restrictions on the boundary operators. 

A chain complex C is defined by a sequence of spaces Cj and by a sequence of maps di, 
which we call boundary operators, where di maps from Cj to Cj_i such that <9j_i<9j = 0. 
Each Ci will be taken to be either Z®* for some given integer p and some given integers Di, 
or more generally we can take Ci = R Di for any ring R. Then, the boundary operator di 
can be written as a matrix with .Dj_i rows and Di columns , with entries of the matrix in 
the given ring. We will be interested in the case that the boundary operator is sparse: each 
row should have 0(1) non-zero entries in it and each column should also have 0(1) non-zero 
entries. Note that in the case of a boundary operator corresponding to a simplicial complex 
or cubical complex or other polyhedral complex, the operator d\ has two non-zero entries in 
every column (since every 1-cell has two 0-cells in its boundary) and the number of non-zero 
entries in a given row is equal to the number of 1-cells attached to a given 0-cells and hence 
is bounded if the local geometry of the complex bounded. Note that if p is prime, Z p = F p 
so the Ci are vector spaces with dimensionality Di. 

It is possible to take products of such hypercomplexes. The product we now describe is 
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known in the math literature as constructing a bicomplex from two chain complexes, and 
then taking the total complex of that chain complex; one important point for us is that 
that process will preserve the sparsity properties of the boundary operators. Given two 
complexes, one complex C given by a sequence Cq,Ci, ... with boundary operators 8$, and 
the other complex C given by a sequence C' ,C[, ... with boundary operators d[, we define 
the product of the two complexes C x C by the sequence 

C ® C , (C <8> C[\ © (Ci ® C ) , - (50) 

That is, the j-th space in C x C is given by 

; „(o c, ,)• (si) 

Let dj be the boundary operators for C x C. We define these boundary operator dj by 
describing its action on vectors in each of the vector spaces C, ® Cj_j for z = 0, j. If f j is 
a vector in and v\ ,• is a vector in C' , we define 

9j (vi ® = (a^) ® + (-i)S is> (52) 

The signs (— l) 1 are chosen so that 

di-A = 0. (53) 

By using this definition for the product of two complexes, one can define the product of 
three or more complexes. The product is associative so that (C x C) x C" = C x (C x C", 
so we simply write products of complexes without parentheses. 



B. Z p Quantum Codes Based on Hypercomplexes 

One important aspect of the above definition of hypercomplexes is that Eq. (153]) holds 
with arbitrary coefficients, including real number coefficients, and not just for Z 2 coefficients. 
If in Eq. (1521) we had not included the sign (—1)*^ then the identity Eq. (|53|) would have 
held with Z 2 coefficients but not otherwise. This will allow us to define Z p codes based on 
hypercomplexes for arbitrary integers p. 

In this section, we will generalize the construction of Ref. 1 in two ways. First, as men- 
tioned, we will extend this construction to Z p codes using our definition of hypercomplexes. 
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Second, we will generalize to hypercomplexes which are the product of three or more hyper- 
graphs. 

Rather than review the construction of Ref. Q, we re-express the construction in other 
language. For the case of a hypercomplex based on the product of two hypergraphs using 
a Z p code with p = 2, with the integer q below equal to 1, this repeats the construction 
of Ref. Q. For other complexes, for n 2, or for q ^ 1, this is a slightly more general 
construction. For example, one could take the product of four hypergraphs to obtain a 
chain complex with spaces Co, C4 and pick q = 2 to obtain a different construction. 

Pick an integer q. We define a code with a number of degrees of freedom equal to 
D q . For each such degree of freedom we have an p-dimensional Hilbert space. We define 
matrices U and V as the "clock" and "shift" matrices. U is a diagonal matrix with entries 
1, exp(2iri/p), exp(47ri/p), ... and V is a permutation matrix with a i in the i,j position of 
the matrix if % = (j + 1) mod p. We define the Hamiltonian by 

H = —A — B, (54) 

where 

A = - A, (55) 
fc=i 

and 



d, 



9+1 



B = - B k - (56) 



fc=i 



In the case that q = 1, the sums of Eqs. fl55|56l) are analogous to the sum over 0-cells and 
2-cells used in Eqs. ( fl5lfT6l) . For example, in Eq. ( 1551) . we sum over k = 1...Dq for q = 1, and 
so we have one term in the sum for each 0-cell. 

Let us introduce some notation. We write matrix elements of the boundary operator <9j 
by {di)i,m, indicating the entry in the Z-th row and m-th column, where we labels the rows 
by integers 1, A-i and the columns by integers 1, Dj. We define A k by 



A k = l[u^ )k ' l + h.c.. (57) 



B k = Hv e {dq+lkk + h.c. (58) 



I 

We define B k by 

-t<>,, i 1 !/.; 

Because <9j_i<9j = 0, we have [A k , Bj\ = for all k, I, for all chain complexes. 
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C. Triviality Properties of Different Hypergraph Codes 



The particular Hamiltonian we considered in Section Ion G x G corresponds to a par- 
ticular example of codes formed in this manner. This code does not have a duality between 
the electric and magnetic degrees of freedom (or rather, the duality that interchanges the 
electric and magnetic degrees of freedom does not have leave the complex unchanged). As 
that imbalance between electric and magnetic degrees of freedom might be responsible for 
the one-sided nature of our result, it is natural to consider instead a more balanced case. 

Suppose we consider a hypergraph which has the same number of hyperedges and hyper- 
vertices. In our language, this hypergraph correspnds to a bipartite graph with the same 
number of vertices in each part, or equivalently to a chain complex where the boundary 
operator d\ is a square matrix. Suppose we choose such a hypergraph at random, with each 
hyperedges containing d hypervertices and each hypervertex in d hyperedges, for some small 
d such as d — 4. Let us fix p to be a prime to make the language simpler: it will enable us 
to talk about the dimension of various vector spaces (fixing p prime is not essential). For 
typical random such hypergraph, the zeroth and first Betti numbers of the complex will be 
0(1). As a result, if we define a complex by taking the product of this complex with itself, 
then the resulting complex will having zeroth, first, and second Betti numbers all being 
0(1). 

The second Betti number counts the "redundancy" of the operators B/,. That is, for any 
v E C2 such that d^v = 0, the operator 

U BV k k ( 59 ) 

k 

is exactly equal to the identity operator; here we define to be the fc-th entry of vector v. 
The second Betti number is the dimension of the space of such vectors because we O3 = 0. 
Similarly, for any vector v in Co such that d[v = 0, 

U A l k ( 6 °) 

k 

equals the identity. The zeroth Betti number counts the "redundancy" of the operators A}.. 
Note that if the complex C had been obtained from a simplicial complex then the zeroth 
Betti number is always non-zero and the product of all the A\. is equal to the identity in 
that case. However, for an arbitrary chain complex C the zeroth Betti number can indeed 
be zero. 
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We now follow an argument due to Yoshidar^ (Yoshida applied this argument to the case 
of Haah's cubic code 1 -). This argument shows the absence of a phase transition as a function 
of temperature for a Hamiltonian with subextensive redundnacy of constraints. Here we give 
this argument for a Hamiltonian with q — 1 given a sequence of complexes for which the 
zeroth and second Betti numbers are both subextensive and for which C3 = 0. 

We first consider the case in which both Betti numbers are zero so that there is no 
redundancy: that is, there is nontrivial product of the operators Ak,Bk that is equal equal 
to the identity. Consider first the case p = 2, where the groups that we use are more familiar. 
Then, there is a qubit for each e and U e is the spin operator S% and V e = Then, we 
can find a unitary Ucuff in the Clifford group to convert the Hamiltonian into the following 
canonical form: 

UcH ff A k Ul, Hff = 2SI (61) 
UcuffB k U cliff = 2Sl +Do . 

Note that the integer k labelling the different A k ranges from 1, D , while the label e in 
S* ranges from l,...,Di. Thus, we map the Dq different operators A k into the operators 
2S% for k = 1,...,Dq and we map the D 2 different operators B k into the operators 2S% +Dq , 
wth k + Dq ranging from Dq + 1, ...,-Do + D 2 . Note that given that both Betti numbers 
are zero, we have Di > D + D 2 . The factor of 2 appears on the right-hand side of the 
above equation because we have defined each operator B^ to equal a product of Pauli 
operators plus its hermitiian conjugate, hence each A kl B k is equal to two times a product 
of Pauli operators. The existence of such a unitary Ucuff follows from the absence of any 
redundancy of the operators, and can be constructed using a process similar to Gaussian 
elimination as follows. Define a rectangular matrix with entries in F 2 , with Dq rows and D\ 
columns, with a 1 in a given matrix element k, I if the operator A^ for k corresponding to 
the given row contains the operator Sf for / corresponding to the given column. For any i, j 
we can find an unitary W in the Clifford group such that WSfW^ = S^S* and such that 
WS^W^ = S% for i k. (This unitary W is simply a controlled-NOT gate in the x-basis). 
So, conjugating with this matrix W changes the entries of the matrix above by adding one 
column to another. Then, by a process of Gaussian elimination, we can use this column 
addition to bring the matrix into the form that matrix element in position k,l is equal 
to 5k t i- Here we use the assumption of no redundancies to show that the original matrix 
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is non-singular. The unitary corresponding to this sequence of row additions changes the 
operators Bk in some way; since the Ak,Bk commute, the operators B k , after conjugation 
by this unitary, are supported only on sites D + 1, ...Di. Then, we can apply a Hadamard 
transformation to make the B k diagonal in the x-basis and then apply a similar Gaussian 
elimination to bring them into the canonical form too. 

For p > 2, the operators U, V generate an extra-special group, rather than the Pauli 
group. The group of normalizers of this extra-special group in the unitary group can be 
regarded as a generalization of the Clifford groups. We can find a unitary in this group of 
normalizers to make a similar transformation to the equaton above, bringing the operators 
into the canonical form: 

U norm A k Ul orm = U k + h.C, (62) 

U norm B k U norrn = Uk+Do + h- c - 
Define the partition function to be 

Z(£) = tr(exp(-/3#)Y (63) 
Define the thermal expectation value of the energy to be 

E((3) = Z(/3)-Ht (h exp(-/3H)) . (64) 



Note that Z(/3) and E(/3) are both invariant under replacing H — » U norm HUl orm , since 
U n0 rm is unitary. However, after the operators A^, Bk are brought into the given canonical 
form, the calculation of the partition function and energy becomes trivial, since everything 
manifestly decouples: for different k, the operators £/& have different support, and we have to 
solve D + D 2 different problems, each problem on a single site. So we can exactly calculate 
E(f3) for this problem. Let 



tr(exp(-/M 

and let 



ti( A k exp(- /3 A k 
E£(P) = ^7 ^, (65) 



E^P) = ^7 ^, (66) 



tr ^exp(— (3Bk 
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Then, 

Z^)-hx[A k e W {-m) = E£{P) (67) 

and 

Z(f3)-hi(B k eM-PH)) = Ei(P). (68) 

That is, the thermal expectation value of A k is the same whether or not the other terms A\ 
for I 7^ k and B\ are present in the Hamiltonian, and similarly for the thermal expectation 
value of .Bfc. As a result, 

Dq-l Dq + 1 

k=l k=l 

The functions E£, E^ can be calculated easily. In fact they are all equal to each other, 
but that fact is not important. In the case p = 2, they are simply hyperbolic tangent 
functions. The crucial thing is that these functions E£(j3), Ej? are all analytic functions 
and are independent of the particular choice of the complex, and so if we consider a family 
of complexes then the thermal expectation value of the energy density 



+ D q+1 



(70) 



is an analytic function and so there is no thermal phase transition as a function of 0. 

Even if the zeroth and second Betti numbers are non-zero, if we consider a sequence 
of complexes in which they are subextensive the limiting thermal expectation value of the 
energy density will still be an analytic function. To see this, remove a subextensive number 
of terms from the Hamiltonian to obtain a new complex with vanishing zeroth and second 
Betti number, and that complex will have an analytic limiting thermal expectation value of 
the energy density. Removing a subextensive number of number of terms will change the 
energy density by o(l) and so will not change the limiting value. 

However, so far as we know, the fact that there is no phase transition in the thermody- 
namics does not necessarily imply that there is no /3 C < oo such that for /3 > (3 C the thermal 



state Z(/3) _1 exp(—(3H) is nontrivial in the sense of Ref. [l(J. This is an interesting problem 
for the future. 

Still, this Hamiltonian is an interesting one to consider. If indeed there is such a (3 C above 
which there is a transition to a circuit nontrivial state, then perhaps this Hamiltonian might 
also have NLTS. Can we find a better candidate Hamiltonian? The Hamiltonian of section 
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H] had a lot of redundancy, which we used crucially in the one-sided proof, but it suffered 
from the fact that the degrees of freedom were on 1-cells and each 1-cell was attached to 
only 2 0-cells. So, another thought is to replace the unbalanced G x G case of section H] with 
the product of a hypergraph with itself, choosing a hypergraph with extensive first Betti 
number so that the resulting complex will still have lots of redundancy (extensive second 
Betti number) and such that the resulting complex will have each 1-cell attached to more 
than 2 0-cells. 

To do even better, a natural guess is that we want a lot of redundancy in both the Ak 
operators and the B^ operators. By choosing an extensive second Betti number, we obtain 
a large redundancy in the operators. Can we obtain an extensive redundancy in both 
and Bk operators? This seems to require going to q > 1. In the discussion later we consider 
whether this is possible; see also Open Question 12.121 

However, homology or cohomology is not the only way to obtain a redundancy in the 
operators. Consider a <i-fold product of a hypergraph chosen to have extensive first Betti 
number. This fc-fold product has extensive ci-th Betti number. Take q — k — 1. Then, there 
is an extensive redundancy in the operators B^. The d — 2-th Betti number may not be 
extensive. However, there can still be redundancy in the operators Ak if d > 2. That is, 
there may still be v G Cd~i such that v is a cocycle. Indeed, choosing v to be any coboundary 
suffices. Thus, a final interesting Hamiltonian to consider is such a <i-fold product, perhaps 
for large d. 

7. DISCUSSION 

Perhaps the most obvious question we have left open is this: can we construct some 
families of complexes and Hamiltonians with NLTS? The Hamiltonians we considered in 
this paper were toric code Hamiltonians. One route to constructing the desired family of 
complexes and Hamiltonians might be to consider other types of Hamiltonians, but let us 
first consider whether we could achieve this goal using toric code Hamiltonians. Let us give 
some intuition on this. We think of expander graphs as being complexes on which particles 
can move very easily: a particle is a O-dimension object and a path is a 1-dimensional object 
and even deleting a sufficiently small nonzero fraction of vertices from an expander graphs 
gives a graph which has long paths. We think of a non-l-hyperfinite complexes as being a 
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complex on a which a string, a 1-dimensional object, can move very easily, as the path of 
a string is a 2- dimensional object and we can delete a small fraction of cells from a non-1- 
hyperfinite complex and get a complex which has large surfaces. Theorem 15.41 considers a 
toric code Hamiltonian on a 2-complex that was not a triangulation of a manifold. However, 
we believe that a similar result could be proven using a (1, 3)-code on the 4-complex which 
is the product of a hyperbolic surface with itself, a complex which is a triangulation of a 
manifold and which is non-1-hyperfmite. Such a (1, 3)-code has degrees of freedom on 1-cells, 
and has the same Hamiltonian as H tc above: there are terms in the Hamiltonian A s for each 
0-cell and B p for each 2-cell and the result that we conjecture would exist in this case is 
that there is no low energy trivial state in which all the A s terms are in their ground state 
(no low energy trivial state without electric defects). Now, suppose we want to come up 
with a (q, r)-code on a ci-complex for which there is NLTS without imposing any restriction 
that there be no electric defects. A natural guess is that we need the complex to be both 
non-g-hyperfmite and non-r-hyperfmite. Perhaps we want the complex to be a triangulation 
of a d- dimensional hyperbolic manifold. See Open Question 12. 121 for the question of whether 
such a complex exists. If d is even, the most favorable situation to achieve this would be 
q = r = d/2 while if d is odd the most favorable situation is q = \d/2\,r = \d/2] or 
q = \d/2],r = \_d/2\. So, we want a <i-complex that is non-[ci/2"|-hyperfinite. 

One motivation to consider complexes which are both non-g-hyperfmite and non-r- 
hyperfmite is based on the result in the appendix of Ref. llCj. There, it was shown that 



for a (1,2) code on a 3-dimensional cubic lattice with a small density of holes removed from 
the lattice there was a trivial ground state. This lattice with a small density of holes removed 
is 2-localizable with a range dependent upon the density of holes This result was based on 
specific properties of this (1,2) code while for more general Hamiltonians it may not be 
case that such a trivial state exists. The result uses the construction of a specific unitary 
quantum circuit for this code. It seems that the ideas in this construction could perhaps be 
extended to show that if the complex is g-localizable or r-localizable that a (q, r) code would 
have a trivial ground state and hence that if the complex is g-hyperfinite or r-hyperfinite 
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that a (q, r) code would have a trivial low energy state. 
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Appendix A: Notions for infinite graphs and noncompact spaces 

Elek— defined a family of graphs {Aj} to be hyperfinite ifVe>03r>0 such that by 
deleting a fixed fraction of edges the resulting components each had diameter < r. 

We have renamed Elek's notion "0-hyperfinite" and made the 
Definition: A family of simplicial complexes (with metrics) {X{\ is /c-hyperfmite if V e > 
3 r > such that by deleting an e-fraction of the simplicies of each X{ to get {Yi}, then 
there are continuous functions /$ : 1$ — > Ki, a /c-complex so that diam(/ i ~ 1 (p)) < r for 
all p G Ki. 

In (arXiv: 1105.320) Abert and Elek showed that the group graphs C of an amenable 
group G (and the Cayley graphs of relatively amenable subgroups) exhibited the hyperfinite 
property "all at once": for all e > 0, there is an r > and an edge set E t of density < e 
such that the components of C\E £ all have diameter < r. 

There is a logical problem with trying to establish the converse: we do not know how 
to even define the e-density of a subset if we do not first know that the graph is amenable. 
(This is representative of a host of issues surrounding amenability, e.g. even in the simplest 
case G = Z constructing the invariant measure requires the axiom of choice.) 

We would like to propose a definition of a concept fc-amenable which can be applied to 
a (single) noncompact metric space X in the favorable case that they possess a transitive 
group of isometries. This definition switches contexts: instead of discussing subsets Y G X 
whose fractional volume vol Y/vol X may make no sense, we discuss invariant measures on 



48 



such subsets. This approach was used by L. Bowen and C. Radin 4 apparently following a 
suggestion of O. Schramm to create a context for studying (radius = r)-sphere packings 
in hyperbolic space. Instead of asking: "which packing is most dense," they ask "which 
translation invariant measure on sphere packings assigns the highest probability to an ar- 
bitrary base point * lying within one of the spheres". We conclude the appendix with our 
definition of fc-amenable, intended as the non-compact or infinite measure counterpart to 
fc-hyperfinite. Just as j-hyperfinite implies fc-hyperfinite for j < k, j-amenable implies k- 
amenable for j < k, and O-amenable should coincide with amenable in cases when both can 
be defined. 

Definition: Let X be a metric space with a transitive group G of isometries. X is k-amenable 
ifVe>03r>0 and a G-invariant u-additive (Borel) probability measure /i on closed 
subsets of X so that: 

1. the probability of a base point iGl lying in a \x random subset is > (1 — e), and 

2. \i is supported on closed subsets admitting continuous maps fa : C{ — > Ki, Ki a 
/c-complex, with diam(/ i _1 (p)) < r for all points p G K^. 
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